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ABSTRACT 



Because of the development of large- for mat, wide-field cameras, microlensing surveys 
are now able to monitor millions of stars with sufficient cadence to detect planets. These 
new discoveries will span the full range of significance levels including planetary signals 
too small to be distinguished from the noise. At present, we do not understand where 
the threshold is for detecting planets. MOA-2011-BLG-293Lb is the first planet to be 
published from the new surveys, and it also has substantial followup observations. This 
planet is robustly detected in survey-l- followup data (Ax 2 ~ 5400). The planet/host 
mass ratio is q = 5.3 ±0.2 x 10 -3 . The best fit projected separation is s = 0.548 ±0.005 
Einstein radii. However, due to the s -H- s^ 1 degeneracy, projected separations of s" 1 
are only marginally disfavored at A% 2 = 3. A Bayesian estimate of the host mass gives 
Ml = 0A3~^q'1y Mq, with a sharp upper limit of Ml < 1.2 Mq from upper limits on 
the lens flux. Hence, the planet mass is m p = 2.4l | g Mj up , and the physical projected 
separation is either rj_ ~ 1.0 AU or r± ~ 3.4 AU. We show that survey data alone 
predict this solution and are able to characterize the planet, but the A% 2 is much 
smaller (A% 2 ~ 500) than with the followup data. The Ax 2 for the survey data alone 
is smaller than for any other securely detected planet. This event suggests a means to 
probe the detection threshold, by analyzing a large sample of events like MOA-2011- 
BLG-293, which have both followup data and high cadence survey data, to provide a 
guide for the interpretation of pure survey microlensing data. 



1. Introduction 

Large- format, wide-field cameras have placed microlensi ng on the cusp of j pining RV and 



transits as a technique able to find dozens of planets at a time (jShvartzvald fc Maozll201ll ). moving 
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the field from the discovery of individual objects to the study of planet populations. Using these 
new cameras, "second generation" microlensing surveys will be able to effectively monitor an order 
of magnitude more events for anomalies due to planets. At the same time, they can maintain an 
observing strategy that makes no reference to whether or not a planetary signal is suspected, thus 
enabling a statistically robust sample of events whose detection efficiencies are well understood. One 
requirement for such a sample is that all events must be analyzed for planets, including signals at 
the limits of detectability. At present the detection threshold is poorly understood since the current 
practice is only to analyze the most obvious signals. However, it is known that microlensing data 
have systematics and correlated noise that make it difficult to use standard statistical measures to 
set the detection thresholds. In this paper, we analyze the microlensing event MOA-2011-BLG-293, 
which is covered by all three second-generation survey telescopes and also has substantial followup 
data, and we suggest a means to study the boundary of what is detectable in the second generation 
surveys. 

Originally, the purpose of microlensing surveys was simply to identify ongoing microlensing 
events, which requires monitoring several million stars with a cadence of about once per night. 
Because a typical planetary signal lasts for only a few hours, it was nearly impossible to detect 
planets from the early survey data. Thus, in order to detect planets, higher cadence followup data 
were needed^. One followup strategy is to monitor one or more targets with increased cadence 
to provide additional coverage of the light curve and to search for anomalies. A second strategy 
is continuous or near-continuous monitoring of a single event of interest, usually because it is 
suspected to be high magnification (A mSuX > 100) or anomalous. These additional observations can 
be taken either by dedicated followup groups or the surveys themselves can go into followup mode 
(typically, continuous or near-continuous observations) if they deem an event to be of sufficient 
interest. Therefore, in followup mode, both survey and followup groups may modify their target 
list and/or observing cadence in response to suspected planetary signals. This strategy has been 
effective at finding planets but makes understanding the dete ction efficienc i es com plex, although 



this has been don e successfully for high magnification events in lGould et al.l (|2010al ). Additionally 



Sumi et al.l (]2010l ) were able to derive a slope (but not the normalization) for the mass ratio function 
of planets from the planetary events known at the time. Of the 13 microl ensing planets publi shed 
to data^l, only one was published from data taken in a pure survey mode ( Bennett et al. 2008l )R 



The new high-cadence, systematic surveys will have sufficient cadence and data quality to de- 
tect and characterize planets with masses as small as the Earth without additional followup data 



Implicit in this is that the data quality is good enough for planet detection. 

''For completeness, we note that there are at least two other planets claimed in the literature, which are not 
considered secure de tections. A possible circumbinary planet proposed to explain anomalies in MAC HO 97-BLG-41 
( Bennett et~ai1ll999h can also be explained by orbital motion of the binary lens (jAlbrow et alj|20pcp. There is also 
evidence for a planetary companion to the lens in MACHO 9 
contains a discussion of why this is inadequate for detection). 



-BLG-35, but only with A X = 20 (|Gaudi et al 



2002 



At least one other pure survey detection should be published soon fsee lBennett et al. 



m pre 



P). 
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(|Gaudill2008l ). Such pure-survey detections require near-24-hour monitoring with a cadence of sev- 
eral observations per hour. Many of these future discoveries will be part of a rigorous experiment 
wherein the detection efficiencies are well understood because they will be found in blind or blinded 
(in which followup data are removed) searches. High cadence surveys, even without g lobal cover 



age, also allow additional science such as the detection of very short timescale events (jSumi et al 



2011 



Recent upgrades by the Optical Gravitational Lensing Experiment (OGLE; Chile) and Mi- 
crolensing Observations in Astrophysics (MOA; New Zealand) collaborations augmented by the 
Wise Observatory (Israel) survey now allow near-continuous monitoring (observations ey ery 15-30 



minutes) of several fields in the Galactic Bulge, 22 hours/day (jShvartzvald &; Maoal2011 



The power of these new surveys comes from the combination of high-cadence, systematic ob- 
servations, which were previously only achievable through followup for a small subset of events, 
and the ability to monitor millions of stars. At the same time, followup observations maintain 
some advantages over current surveys. Because the followup networks have access to additional 
telescopes at various sites, followup observations often have redundancy. This makes them less 
vulnerable to bad weather, which can create gaps in the data. Additionally, multiple data sets 
at a given epoch provide a check on s ystematics or other as trophysical phenomena that may cre- 
ate false microlensing-like signals (see iGould et al. Ilin prepl ). Simultaneous or near simultaneous 
observations from multiple sites are also required to measure terrestrial microlens parallax (e.g. 
Gould et al. 1 12003 ) . Furthermore, since followup observations are targeted, they can achieve a 
much higher cadence, and are frequently continuous, although the current strategy for surveys is 
typically to switch to near continuous followup observations for events of interest. Finally, while 
followup groups routinely make an intensive effort to get observations in additional filterqj, survey 
groups are less aggressive about obtaining such observations. All of these additional bits of informa- 
tion can increase confidence in the microlensing interpretation and reduce ambiguity in the models. 
The trade off is that with somewhat sparser data coverage, surveys are able to systematically 
monitor more than an order of magnitude more events. 

The additional planets detected by surveys, which are not currently being detected with fol- 
lowup, will fall into two categories related to the two kinds of caustics produced by a 2-body lens. 
First, there will be many more planetary caustic anomalies detected. These caustics are created 
along or near the planet-star axis at a distance from the lens star that depends on the projected 
separation, s, and the mass ratio, q, between the two bodies. Anomalies created by these caustics 
can be found with current followup but since the source trajectory is random with respect to the 
binary axis, these anomalies occur in a random place in the light curve. Hence, surveys will detect 
more of them because they can observe more stars. Second, a planetary companion to the lens 
induces a caustic at the position of the lens star, the so-called "central" caustic. Source crossings 
of such caustics can be predicted in advance because they require that the source trajectory pass 



5 Microlensing observations are normally done in /-band (or similar filters) because that is the optical band that 
is most sensitive toward the Galactic Bulge. In order to derive source colors as in Sec. [3J we need observations in a 
different filter, for which we typically use V-band. 
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very close to the position of the lens star. Surveys will observe more central caustic events that are 
too faint to observe with current followup or are not recognized to be high magnification quickly 
enough to organize followup observations. For both types of events, there is the question of whether 
the survey data alone are indeed sufficient to detect planets in individual microlensing light curves 
in spite of having sparser data. For central caustic events there is an additional question of whether 
or not the anomaly will be sufficiently characterized, since the models can be quite degenerate for 
these kinds of events, sometimes with little constraint on the mass ratio between the lens star and 



its companion (e.g. IChoi et al.ll2012l ). Given this much larger sample of events which will contain 
signals of all significance levels including ones that can be confused with systematics, the challenge 
is to create a subset of events for which the vast majority of planetary signals can be considered 
reliable, and secondarily for which the planets are well-characterized, in the case of central caustic- 
type events. I Gould et al.1 (j2010al ) estimate a threshold of A% 2 =350-700 would be appropriate, 



but the true value is unknown. In principle, such questions could be addressed with simulations. 
However, in simulations it is difficult to account for real effects such as data systematics and stellar 
variability. Hence, using actual microlens data provides field testing that complements results from 
simulations. 

MOA-2011-BLG-293 provides an opportunity for investigating survey-only detection thresh- 
olds. The planet is robustly detected in the survey-l- followup data (A% 2 ~ 5400), and the event was 
observed by all three current survey telescopes. Wise Observatory obtained data of the anomaly 
in their normal survey mode without changing their observing cadence, and the rest of the light 
curve is reasonably well covered by OGLE and MOA survey data. For this event, we are able to 
determine whether the survey data alone can successfully "predict" the solution determined when 
all of the data are included. 

This event also has the faintest source of any published planetary event. We show that for such 
faint sources small systematic errors in the flux measurements can radically affect the microlensing 
solution, even when all the anomalous features occur at high magnification when the source is 
bright. In particular, the source flux and the event timescale are determined primarily from data 
near baseline where small systematic errors may be of order the change in flux being measured. 
Because the systematic errors in the timescale propagate to many other quantities including the 
planet/star mass ratio, they must be investigated carefully. This is particularly important for future 
surveys where many of the events will be at or beyond the magnitude limit at baseline. 

We begin by presenting the discovery and observations of MOA-2011-BLG-293 in Section [2j 
The color-magnitude diagram of the event is presented in Section [3] and used to derive the intrinsic 
source flux. In Section [U we address the consequences of systematics in the measured flux when 
they are similar in magnitude to the source flux. Then, in Section [5] we present the analysis of the 
light curve of the event, and we compare the results with and without followup data in Section 
15.21 Additional properties of the event are derived in Section [6l and the physical properties of the 
lens star and planet are derived from a Galactic model in Section [7J We discuss the implications 
for future survey-only detections in Section [BJ Finally, the possibility of detecting the lens with 
adaptive optics (AO) observations is discussed in the Appendix. 
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2. Data Collection and Reduction 



MOA issued an electronic alert for MOA-2011-BLG-293 [(RA,Dec) = (17:55:39.35, -28:28:36.65), 
(Lb)=(1.52 -1.66)] at UT 10:27, 4 Jul 2011 (HJD' = HJD-2450000 = 5746.94), based on survey ob- 
servations from their 1.8m telescope with a broad R/I filter and 2.2 deg 2 imager at Mt. John, New 
Zealand. At UT 12:45, the Microlensing Follow-Up Network (/iFUN) refit the data and announced 
that this was a possible high-magnification event, where "high-magnification" is A m3iX > 100. At 
UT 17:28, //FUN upgraded to a full high -magnification alert (-^max 

> 270), emailing subscribers to 

their email alert service, which includes members of /xFUN and other microlensing groups, to urge 
observations from Africa, South America, and Israel. Additionally, a shortened version of the alert 
was posted to Twitter. This prompted /xFUN Weizmann to initiate the first followup observations 
at UT 19:45, using their 0.4m telescope (/ band) at the Martin S. Kraar Observatory located on 
top of the accelerator tower at the Weizmann Institute of Science Campus in Rehovot, Israel. At 
UT 23:25, /xFUN Chile initiated continuous observations using the SMARTS 1.3m telescope at 
CTIO. At UT 00:00 /xFUN issued an anomaly alert based on the first four photometry points from 
CTIO, which were rapidly declining when the expected behavior was rapid brightening. The great 
majority of the CTIO observations were in I band, but seven observations were taken in V band 
to measure the source color. In addition, the SMARTS ANDICAM camera takes H band images 
simultaneously with each V and / image. These are not used in the light curve analysis but are 
important in the Appendix. 

MOA-2011-BLG-293 lies within the survey footprint of the MOA, OGLE, and Wise microlens- 
ing surveys and so was scheduled for "automatic observations" at high cadence at all three observa- 
tories. MOA observed this event at least 5 times per hour. Wise observed this field 10 times during 
the 4.6 hours that it was visible from their 1.0m telescope, equipped with 1 deg 2 imager and /-band 
filter, at Mitzpe Ramon, Israel. The event lies in OGLE field 504, one of three very high cadence 
fields, which OGLE would normally observe about 3 times per hour. In fact, it was observed at 
a much higher rate, but with the same exposure time, in response to the high-magnification alert 
and anomaly alert. Unfortunately, high winds prevented opening of the telescope until UT 01:02. 
OGLE employs the 1.3m Warsaw telescope at Las Campanas Observatory in Chile, equipped with 
a 1.4 deg 2 imager primarily using an /-band filter. 

The data are shown in Figure HJ Several features should be noted. First, the pronounced part 
of the anomaly lasts just 4 hours beginning at HJD' = 5747.40. The main feature is quite striking, 
becoming about one magnitude brighter in about one hour. The coverage during the anomaly is 
temporally disjoint between the observatories in Israel and those in Chile, a point to which we 
return below. Finally, the CTIO data show a discontinuous change of slope ("break"), which is the 
hallmark of a caustic exit, when the source passes from being partially or fully inside a caustic to 
being fully outside the caustic (see Fig. 0]). 



MOA and OGLE data were reduced using their standard pipelines (jBond et al.ll200ll ; lUdalski 



20031 ) which are based on difference image analysis (DIA). In the case of the OGLE data, the source 



is undetected in the template image. Since the OGLE pipeline reports photometry in magnitudes, 
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an artificial blend star with a flux of 800 units (Iogle = 20.44) was added to the position of the 
event to prevent measurements of negative flux (and undefined magnitudes) at baseline when the 
source is unmagnified. 



Data from the remaining three observatories were also reduced using DIA (jWozniakl |2000j) , 
with each reduction specifically adapted to that imager. Using comparison stars, the Wise and 
Weizmann photometry were ali gned to the same flux scale as the CTIO I band by inverting the 



technique of I Gould et al.l (|2010bl ). That is, the instrumental source color was determined from CTIO 
observations, and then the instrumental flux ratios (CTIO vs. Wise, or CTIO vs. Weizmann) were 
measured for field stars of similar color. The uncertainties in these flux alignments are 0.016 mag 
for Wise and 0.061 mag for Weizmann. 

2.1. Data Binning and Error Normalization 

Since photometry packages typically underestimate the true errors, which have a contribution 
from systematics, we renormalize the error bars on the data, as is done for most microlensing 
events. After finding an initial model, we calculate the cumulative x 2 distribution for each set of 
data sorted by magnification. We renormalize the error bars using the formula 



Vi = k \Jv 2 + <in (!) 

and choosing values of k and e mm such that the x 2 per degree of freedom x 2 cc j = 1 and the cumulative 
sum of x 2 is approximately linear as a function of source magnification. Specifically, we sort the 
data points by magnification, calculate the A% 2 contributed by each point, and plot Yli ^xf 
as a function of N to create the cumulative sum of x 2 > where TV" is the number of points with 
magnification less than or equal to the magnification of point N. Note that <7j is the uncertainty 
in magnitudes (rather than flux). The values of k and e m i n for each data set are given in Table 
[TJ Except for OGLE, the values of e m i n are all zero. This term compensates for unrealistically 
small uncertainties in the measured magnitude, which can happen when the event is bright and the 
Poisson flux errors are small. 

For the MOA data, we eliminate all observations with t outside the interval 5743.5 < t(HJD') < 
5749.5 (see Section l4T2j) . We also exclude all MOA points with seeing > 5" because these data show 
a strong nonlinear trend with seeing at baseline. After making these cuts, we renormalize the data 
as described above. 

To speed computation, the OGLE and MOA data in the wings of the event were binned. In 
the process of the binning, 3a outliers were removed. This binning does not account for correlations 
in the data, which if they exist can increase the reduced x 2 above the nominal value of x 2 c d = 1- 

3. CMD 

We use the CTIO / and V band data to construct a color-magnitude diagram (CMD) of 
the event. We measure the instrumental (uncalibrated) source color by linear regression of the 
V and / fluxes (which is independent of the model) and the magnitude from the fs,CTio of our 
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best-fit model: (V — I,I)s = (0.37, 22.27) ± (0.03, 0.05). The position of the source relative to the 
field stars within 60" of the source (small dots) is shown in Figure [2] as the solid black dot. We 
calibrate these magnitudes and account for the reddening toward the field by assuming the source 
is in (or at least suffers the same extinction as) the Bulge and using the centroid of the red clump 
as a standard candle. Because of strong differential extinction across the field, we use only stars 
within 60" of the source to measure the centroid of the red clump. Since the event is in a low 
latitude field, there are more stars than is typical for bulge fields and the red-clump centroid can be 
reliably determined even with this restriction. In instrumental magnitudes, the centroid of the red 
clump is (V — I,I)r i = (0.59, 16. 90) compared to its intrinsic value of (V — I, I) c i,o = (1-06, 14.32) 
(jBensbv et al.l l201ll : iNatafl 120120 , which assumes a Galactocentric distan ce Ftp = 8kpc an d that 
the mean clump distance toward I = 1.5 lies 0.1 mag closer than Rq (jRattenburvl 120071 ) . We 
can apply the offset between these two values to the source color and magnitude to obtain the 
calibrated, deredden ed values (V — I,I)g q = (0.84,19.69) ± (0.05,0.16). The uncertainty in the 
color is derived from iBensbv 



by comparing the spectroscopic colors to the microlens 
colors of that sample. The uncertainty in the calibrated magnitude is the sum in quadrature of 
the uncertainty in fs, ctio from the models (0.05 mag), the uncertainty in Rq (5%— >0.1 mag), the 
uncertainty in the intrinsic clump magnitude (0.05 mag), and the uncertainty in centroiding the 
red clump (0.1 mag). 



4. Effect of Faint Sources on Microlens Parameters 



The source star in MOA-2011-BLG-293 is extremely faint with an apparent magnitude in the 
OGLE photometry of is, ogle = 21.7. Consequently, the measured flux errors can be comparable 
to or larger than the source flux, particularly near baseline. Because of this, systematics in the 
baseline data must be carefully accounted for so as not to bias the microlens results. Systematics 
in the measured flux at the level of fs can lead to biases in the measured Einstein timescale, ts, 
of the same order. We begin by discussing robust parameters, which can be measured solely from 
the highly magnified portion of the light curve and so, are independent of uncertainties in the flux 
measured near baseline. Then in Section 14.21 we discuss in detail the effect of systematics in the 
measured baseline flux on the microlens parameters, particularly t& and the mass ratio between 
the components of the lens, q. 



4.1. Robustly Measured Parameters 



At high magnification, a microlensing light curve for a point source being lensed by a point 
lens can be described by the unmagnified, baseline flux of the event, /base) an d three parameters 
("invariants"): the time of the peak, tp . the differen ce between /base and the peak flux, /ii m , and the 
effective width of the light curve, t e g (|Gouldlll996l . see Eq. (2.4) and (2.5) in). These parameters 
are nearly invariant under changes to the source flux and are robustly determined by the light 
curve. The change in the observed flux due to the event can then be written as a function of these 
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invariants: 

/obs(*) — /base = G^(t; t , t c g, fn m ), (2) 

where / Q b s is the observed flux and the subscript on G refers to the number of parameters. Note 
that /base is also an observable. 

In the limit where the event is highly magnified, the exact form of Gs(t) can be derived from 
the microlens variables. The three microlens variables of a point-source-point-lens microlensing 
model are to, the impact parameter in units of the Einstein radius, uq, and the Einstein crossing 
time, £e- The observed flux is given by 

/obs = f S A(t) + f B = f s [A(t) - 1] + / base , (3) 

where A(t) is the magnification of the source, fs is the flux of the source, and fs is the flux of all 
other stars blended into the PSF (including the flux from the lens). By definition, /base = fs + Ib- 
For a point lens in the limit of high magnification (A(t) 3> 1), the magnification is given by 

A ® - — Attt ( 4 ) 

uoQ(t) 

where 



Q(t;to,*eff) = yi+(^) , (5) 

is a function of only time and the invariant: 

*eff = UotE- (6) 
In this limit, the evolution of the observed flux, Gs(t; to, t e g, /n m ), is then given by 

G 3 (t;to,t cS ,f hm ) = where /nm = — . (7) 

Q(t) uo 

If finite source effects are detected, the change in the observed flux is a more complicated 
function because of the additional microlens parameter p, which is the source size in units of the 
Einstein radius. However, there is also an additional invariant 

t* = ptE, (8) 

the source crossing time, which determines the width of the peak of the light curve for a point lens. 
Hence, the change in the observed flux can be written as 

fobs(t) — /base = G±{t\ t ,t o g, fu m , t+) = G${t)B {Qt c g/t*) , (9) 

w here B(Qt P ff/ t+) = B(u/p) is a funct ion composed of elliptic integrals, whose exact form is derived 



m 



Gould (119941 ) and lYoo et al.l ((2004J). 



In the case of a two-body lens like MOA-2011-BLG-293, the invariants may not be obvious 
from the light curve, but they are still robustly measured as we show below. For example, the 
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width of the peak is distorted by two-body perturbation, but based on the source trajectory (Fig. 

it can be seen that the width of the first bump at HJD'5747.46, which is caused by the cusp 
crossing, will be slightly larger than 2i*. For a two-body lens with a central caustic crossing, there 
is also another invariant due to the mass ratio, q, between the two lensing bodies, so 

fobs(t) — /base = G^'i *0, *eff j /lim> **> <7*e)- (10) 

This new invariant qt-^ can be understood as follows. For central caustic^, like the one in MOA- 
2011-BLG-293, the caustic size is proportional to the mass ratio of the two lensing bodies, q, and 
the caustic shape is roughly constant for a given s. Therefore, the time between successive features 
in the light curve is set by qt~E, i.e., the size of the caustic multiplied by the characteristic timescale, 
and since the observed times of the features can be well measured, the uncertainty in this quantity 
is extremely small. In this case, the main features are the two bumps in the light curve and the 
discontinuity in the slope that occurs between the bumps. [Note that a two-body lens introduces 
two parameters in addition to q: the separation between the two bodies projected onto the plane 
of the sky, s, and the angle of the source trajectory with respect to the binary axis, a. Because 
these parameters are of less interest, we do not discuss them in this context.] 

Table [2] shows that for this event these quantities, t e g, /ii m , £*, and qt&, do indeed have 
extremely small uncertainties and can approximately be considered invariants. 



4.2. Parameters Vulnerable to Systematics 

The above invariants are determined by the data taken near the peak of the light curve where 
A S> 1. However, in order to extract the values of the microlens parameters, uq, p, and q, from 
the invariants, we must measure ig. The information onfg must neces sarily come from the wings 



of the light curve where the magnification is small and A(t) ^ l/u(t) (jDominikl 120091 ) . Since the 
magnification is small, the change in the observed flux compared to the source flux is also small, 
so the measurement of tE may be considered to be the statistical sum of many measurements of a 
small change in flux. In order to get an accurate measurement of ts, the statistical and systematic 
errors in the flux must be smaller than the change we are trying to measure. 

In the case of MOA-2011-BLG-293, because the source flux is extremely faint, it is difficult 
to measure accurately. When the magnification is a factor of a few or less (in the wings and at 
baseline), if the flux is not measured with an accuracy substantially smaller than the source flux, 
this can lead to bias in the measurement of fs or equivalently t-p,, since fstE = /lim^eff is robustly 
determined, and so to bias in quantities dependent on ts such as q. To check for this possible 
source of bias, we bin the OGLE and MOA data by 30 days to see if their measurements of the 
baseline flux are stable (Fig. [3]). We find that the OGLE measurement of the baseline flux is stable 



Two-body lenses with unequal mass ratios will create one caustic at the position of the more massive body, the 
"central" caustic, and another set of caustics elsewhere, the "planetary" caustics. 
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at a level that is smaller than the observed source flux. Therefore, we use all of the OGLE data in 
our models. We note that the flux after the event (t > HJD'5790) appears to be at a lower level 
than the baseline before the event. In Section [5.14 we discuss the effect of assuming the baseline 
decreases at a constant rate during the course of the event. 

As seen in Figure [3l the MOA baseline flux exhibits scatter in excess of the measured pho- 
tometric errors, and there is also variation in measured baseline flux from season to season. The 
magnitude of this scatter is similar to the magnitude of the source flux. Because of this vari- 
ation, we conclude that the baseline flux is not sufficiently well measured in the MOA data to 
detect the small changes in flux necessary to measure t^- As a result, to avoid biasing our results, 
we use only the MOA data from the peak of the light curve where the photometry is precise: 
5743.5 < i(HJD') < 5749.5. 

5. Analysis 



Without any modeling, we can make some basic inferences about the relevant microlens param- 
eters from inspection of the light curve. MOA-2011-BLG-293 increases in brightness from / ~ 19.7 
to I ~ 15.0, indicating a source magnification of at least 75. Additionally, except for the deviations 
at the peak, the event is symmetric about to- From these two properties, we infer that only central 
or resonant caustics (both of which are centered on the position of the primary) are relevant to the 
search for microlens models. 

We fit the light curve using a Markov Chain Monte Carlo (MCMC) procedure. In addition to 
the parameters described in Section^ a model with a two-body lens has two additional parameters: 
the angle of the source trajectory with respect to the binary axis defined to be positive in the 
clockwise direction^, a, and the projected separation between the two components of the lens 
scaled to the Einstein radius, s. Because they are approximately constants, we use the parameters 
t e s and t* in place of the microlens variables uq and p. For a given model, Equation ([3]) must 
be evaluated for each observatory, i, so fs, fo — > f,s>;, £b ,i- We adopt the "natural" linear limb- 



darkening coefficients T = 2u/(3 — u) (lAlbrow et al 



source in the CMD, we estimate that T c g 
darkening coefficients for T c q = 5250K and T c g 
s^ 1 to find r y = 0.6368 and Tj = 0.4602. 



19991 ). Based on the measured position of the 
5315K and log g = 4.5 cg s . We average the linear limb- 
5500K from IClaretl (|2000i ) assuming v tur b = 2 km 



The magnifications are calculated on an (s, q) grid, using the "map-making " technique (IDong et al 



2009 



2006 ) in the strong finite-source regime and the "hexadecapole" approximation (jPejcha &: Heyrovsky 



Gouldll2008l ) in the intermediate regime. 



We began by searching a grid of s and q to obtain a basic solution for the light curve. For 
central caustic crossing events like this one, there is a well known degeneracy betwee n models with 
close topologies (s < 1) and wide topologies (s > 1) (e.g. iGriest k, Safizadehlll998l ), We initially 
searched a broad grid for close topologies and then used the results to inform our search for wide 



7 The binary axis has its origin at the center of magnification and is positive in the direction of the planet. 
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solutions, since to first order, s — > s . The basic model from this broad grid has s ~ 0.55, 
q ~ 0.005, and a ~ 220°, such that the source passes over a cusp at the back end of a central 
caustic. This caustic is created by a two-body lens with a mass ratio similar to that of a massive 
Jovian planet orbiting a star. Figure |4] shows this basic geometry with the source trajectory relative 
to the caustic structure. The bump in the light curve at HJD' ~ 5747.45 is created when the source 
passes over the cusp of the caustic. 

Because the Wise and Weizmann data only overlap with other data sets where their errors are 
extremely large, there is some concern that the parameters of the models will be poorly constrained, 
since within the standard modeling approach the flux levels of these data can be arbitrarily adjusted 
up or down relative to the other data. However, from the flux alignment described in Section [2] we 
have an estimate of fsi for these data relative to fs,CTio- This alignment gives us an independent 
means to test the validity of our model. If the model is correct, then the values of fs,wise an d 
7s,Weizmann should agree with /s,ctio within the allowed uncertainties. Alternatively, if we include 
the flux- alignment constraint in the MCMC fits, the solution should not change significantly. 

We incorporate the flux-alignment constraint in a way that is parallel to the model constraints 
from the data, i.e., by introducing a \ 2 penalty: 

2 (/s.CTIO - fs,i) 2 111 10 //s.CTIO + fs,i\ /-.-.x 

Xb = \ <~ ; aam ' i = ^5\ 2 )**' 

where i corresponds to the observatory with the constraint, and <7j is the uncertainty in magnitudes 
of the flux alignment for that observatory. In the absence of any constraints, the flux parameters 
for each observatory, fg t i and fs,i, are linear and their values for a particular model can be found 
by inverting a block-diagonal covariance matrix, b. We include the flux constraints by adding half 
of the second derivatives of x\ to the b matrix: 

A6(/s,i,/5,fc) = ^j—, (12) 

where k =CTIO and 5ik is a Kronecker-delta. This couples formerly independent 2x2 blocks. 
Strictly speaking, the equation for <7fl UXj j given in Equation (|11|) is a numerical approximation. 
Therefore, we iterate the linear fit until the value of <7fl UXj j is converged, which typically occurs in 
only a few iterations. 

We refined the (s,q) grid around our initial close solution, fitting the data both with and 
without flux- alignment constraints. The mean and 1 a confidence intervals for the parameters from 
these two fits are given in Table El There are only small quantitative differences between the two 
solutions, and nothing that changes the qualitative behavior of the model. The slight increase in 
X 2 is expected because of the additional term due to the flux constraints. After finding this close 
solution, we repeated the grid with s — > s _1 to identify the wide solution. The parameters of this 
solution are also given in Table [3] both with and without flux- alignment constraints. The close 
solution is mildly preferred over the wide solution by Ax 2 ~ 3, so we quote the values for the 
flux-constrained close solution: 

q = 5.3 ±0.2 x 10~ 3 s = 0.548 ± 0.005, (13) 
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noting that the two topologies give very similar solutions (except s — > s 1 ). 

Additionally, we searched for a parallax signal in the event by adding two additional free 
parameters to the fit for the close solution: tte,n and tte,e, the North and East components of the 
parallax vector (e.g., iGouldl 12004! ) . The parameters of this fit are given in Table [3l No parallax 
signal was detected, and we found no interesting constraints on these parameters. The \ 2 improves 
for fits including parallax by only A% 2 = 7 for two additional degrees of freedom. In some cases, 
even when parallax is not detected, meaningful upper limits can be placed on the parallax, but in 
this case we have an uninteresting 3a constraint of < |tte| < 7.8. 



5.1. Effect of Systematics in /base 

As discussed in Section f4. 21 it is possible that the underlying OGLE baseline flux is changing 
during the course of the event. In order to test how that could introduce systematic effects in 
our results, we create a fake OGLE data set accounting for a constant decrease in baseline flux 
during the event. Specifically, we assume that the baseline flux decreases at a constant rate between 
HJD'5710. and HJD'5790. leading to an overall decrease in flux of 0.27/5 i ogle- We then repeat 
the MCMC procedure for the close solution including flux- alignment constraints. We find that the 
value of ijj increases by 15%, and consequently, the values of q, uo, and p decrease by the same 
amount. In principle, this could represent a systematic error in our results. However, at the present 
time, the evidence for a change in the baseline flux is weak, so we only report these results for the 
sake of completeness. 

5.2. Analysis with Survey-Only Data 

From this analysis, we have a robustly detected planet (A% 2 ~ 5400 compared to a point lenj§) 
and a well-defined solution. Now we can ask whether the planet could have been detected from the 
survey data alone, whether the solution is well-constrained, and most importantly, whether it is the 
same solution. To begin, we construct a "survey only" subset of the data. We first eliminate the 
Weizmann and CTIO data. Second, we "thin out" the OGLE data to mimic OGLE survey data as 
they would have been if there had been no high-magnification or anomaly alerts. OGLE data on 
several nights previous to (and following) the peak have a cadence of 1 observation per 0.015 days. 
We therefore adopt a subset of 18 (out of 44) OGLE points from the peak night with this sampling 
rate. 

We repeat the analysis on this survey-only data set beginning with a broad grid search and 
then refining the solution following the same procedure used for analyzing the complete data set. 



8 Note that the numbers quoted for the point lens models include constraints from the flux alignment in the 
fit. Removing the flux- alignment constraints improves the x 2 > primarily because the Weizmann data can be scaled 
arbitrarily. However, compared to the planet fit, the point lens fit without flux-alignment constraints is still extremely 
poor, A\ 2 ~ 4400. Flux-alignment constraints have very little effect on the point lens fit to survey only data. 
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We find that even without flux- alignment constraints, the global search isolates solutions in the 
general neighborhood of the solution found from the full data set. The fits to the survey-only data 
set are compared to fits with all data in Figure [5j Here, the Ax 2 of the fit compared to a point lens 
fit for the survey-only data is 487, nearly all of which comes from data in the time-span shown in 
Figure [5] 5747.1 < t(HJD') < 5748.8. This is smaller than the Ax 2 of any published microlensing 
planet. However, the parameters of the fit are well constrained with errors only a factor of 1.5-2 
larger compared to fits with the full data set. Applying the flux-alignment constraint to this model 
confirms its validity, i.e., it does not appreciably change the solution (see Table [3|). In this case, it 
is clear that the survey data are sufficient to robustly detect and characterize the planet. 

In order to push farther into the limits of detectability, we also analyze this event without the 
Wise data, since those data contain most of the deviation from the underlying point lens. With 
only the MOA and thinned OGLE data, we find A% 2 ~ 70 between the 2-body lens and point lens 
models. We note that the point lens model has an unreasonably large value of parallax, -zte ~ 20, 
making it somewhat suspicious. Without parallax, the Ax 2 between the point lens and 2-body fit 
increases to Ax 2 ~ 170. Although this is a factor of 3 smaller than the Ax 2 with the Wise data, 
the constraints on the planet-star mass ratio are still broadly confined to be planetary, assuming 
Ml < O.5M , with a 3a range of 0.001 < q < 0.025. However, because of the small A% 2 , it appears 
unlikely to us that a planetary detection would have been claimed from solely the MOA and OGLE 
data even though the solution is formally well-constrained. 

6. Physical Properties of the Event 

Since finite source effects are measured in this event, we can determine the angular size of the 
Einstein ring, 8~e, and the lens-source relative proper motion, /it. First, we estimate the angular size 
of the source, 6±, from the observed co lor and magnitude. W e transform the (V — I)s,o color to 



(V — K) using the dwarf relation fromlBessell &; Brettl (|1988l ). Then we use the (V — K) surface 



brightness relations from iKervella et al.l (]2004l ) to find 8+ = 0.42 ± 0.03 ^as. From this we derive 



the lens-source relative proper motion and angular Einstein radius, 

H= — = 4.3±0.3masyr~ 1 ; 8 E = /z*e = 0.26 ± 0.02 mas. (14) 

The uncertainties in these quantities come from a variety of factors. Specifically, the uncertainties 
in the Galactocentric distance, Rq, and the measured intrinsic brightness of the red clump, the 
centroiding of the red clump from the CMD, and uncertainty in the surface brightness relations. 
The uncertainty contributed by the surface brightness relations is 0.02 mag, and the uncertainties 
from the other factors are giv en in Section [3 The contribution of these factors can be understood 



from their relationship to 8* (jYee et al.l 12009 ): 



(15) 



where fs is the source flux from the microlensing model and Z captures all other factors. Taking 

's 



account of all factors mentioned above, we find a(Z)/Z = 8%. Since the statistical error in fV 2 



15 



is only 2.3%, the error in Z completely dominates the uncertainty in Q*. In general, the error in 



fs propagates in opposite directions for #e and \i (jYee et al.ll2009l ). However, in the present case, 
since this error is small, the fractional error in these quantities is simply that of Z, as indicated in 
Equation (jl4|) . 

7. Properties of the Lens 

7.1. Limits on the Lens Brightness 

We can use the observed brightness of the event to place constraints on the lens mass. Since 
the source and lens are superposed, any light from the lens should be accounted for by the blend 
flux, fs,i, which sets an upper limit on the light from the lens. The unmagnified source is not 
seen in the OGLE data. From examination of an OGLE image at baseline with good seeing, we 
estimate the upper limit of the blend flux to be Ib,o > 17.77 based on the diffuse background light 
and assuming that the reddening is the same as the red clump. Assuming all of this light is due to 
the lens, the absolute magnitude of the lens is 

M I>L > I B ,o + (A I)S ~ A ItL ) - 5 log = 3.25 + (A I>S - A JjL ) + 5 log (16) 

lUpC JJl 

where Aj^ and Aj^l are the reddening toward the source and lens, respectively, and Dl is the 
distance to the lens. Since the lens must be in front of the source, we have Ai t s > Ar,£- Moreover, 
the lens should be closer than Rq (or at any rate, not much farther). Hence, Ml l > 3.25 is a 



conservative lower limit. From the empirical isochrones of lAn et al.l (|2007l ). this absolute magnitude 
corresponds to an upper limit in the lens mass of Ml < 1.2M . We conclude from these flux- 
alignment constraints that either the lens is a main sequence star or, if it is more massive than our 
upper limit of 1.2M , then it must be a stellar remnant such as a very massive white dwarf or a 
neutron star. 

We can use our measurement of #e to estimate the distance to the lens based on its mass: 

Di - ( S~ L w + ' withKS as? - "■"—"s 1 - < 17 > 

where Ds is the distance to the source. If we assume the source is at 8 kpc (i.e., about 0.1 mag 
behind the mean distance to the clump at this location) and Ml, = 1.2Mq, we find Dl = 7.6 kpc. 
Hence, the lens could be an F/G dwarf or stellar remnant in the Bulge, or it could be a late-type 
star closer to the Sun. 



7.2. Bayesian Analysis 



Similar to lAlcock et al.l (119971 ) and iDominikl (120061 ). we estimate the mass of the lens star 
and its distance using Bayesian analysis accounting for the measured microlensing parameters, the 
brightness constraints on the lens, and a model for the Galaxy. The mathematics are similar to what 
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is described in Section 5 of iBatista et al.l (|201ll ). although the implementation is fundamentally 
different because we do not have meaningful parallax information. Specifically, we perform a 
numerical integral instead of applying the Bayesian analysis to the results of the MCMC procedure. 
We begin with the rate equation for lensing events: 



d 4 F 



dD L dM L d 2 fi 



u(x, y, z)(2R E )v re if(fi)g(M L 



(18) 



where v(x, y, z) is the density of lenses, -Re is the physical Einstein radius, v IC \ is the lens-source 
relative velocity, /(//) is the weighting for the lens-source relative proper motion, and g(Mi) is 
the mass function. The vector form of the lens-source relative proper motion is n, which can be 
described by a mag nitude, /u, and an angle, cj), such that d 2 /j> = iid^d<j). We transform variables (see 
Batista et al.ll201ll ) to find 



d 4 T 



dD LdO^dtftdcj) 



At7T re l 



v(x,y,z)f((j,)g(M L ). 



(19) 



To find the probability density functions for the lens, we integrate this equation over the variables 
9e and 0, using a Gaussian prior for Be with the values given in Eq. (|14p and a flat prior for <f>. 
We calculate /i from t E and #e using Equation (|14p . We also integrate over D$, which appears 
implicitly in 7r re i and /(//). For D$, we include a prior for the density of sources based on our 
Galactic model (see below) assuming the source is in the Bulge. 



we as sume g(M) oc M 
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(20) 



Note that the variables in /„ are given in Galactic coordinates rather than Equatorial coordinates. 
The transformation between the two is simply a rotation by 60°. Still working in Galactic coordi- 
nates, the expected proper motion, /x exp , takes into account the typical motion of a star in the Disk, 
v, and the motion of the Earth during the event, v e = (v® t N sal ,v®,E sal ) — (—0.80, 28.52) km s _1 , 



V L - ( v + v ffi) v 5 - (V© + Vg 



cxp 



D L 



D S 



(21) 



where v 
lOkms" 1 ) 



= (7, 12) km s + (0, -y ro t) and v Iot = 230 km s . For the Disk we use v = (0, v rot — 
and <x = (^,^1, o"/i,£;g al ) = (20, 30) kms -1 , and for the Bulge v = (0, 0)kms _1 and 
f = (cM.-NgaP^M.-EgJ = ( 100 ' 100) kms" 1 . 



For the stellar density v(x,y,z), we use the model from lHan Sz Gouldl (|2003l ) including a bar 
in the Bulge. We assume the Disk has cylindrical symmetry with a hole of radius 1 kpc centered 



'We will neglect constants of proportionality as they are not relevant to a likelihood analysis. 
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at Rq = 8 kpc. We limit the Bulge to 5 < D < 10 kpc, where D is the distance from the observer 
along the line of sight. 

For the Bayesian analysis, we use tE = 21.7 days measured from the microlensing fit to the 
light curve. We also have the constraint from the lens brightness that Ml = 0g/(nTr vc \) < 1.2Mq. 
This analysis implicitly assumes that the lens is a main sequence star. The lens could be a stellar 
remnant, although this is much less likely because of their smaller relative space density. The 
possibility that the lens is a stellar remnant could be tested several years from now when the 
source and lens have moved sufficiently far apart so as to be separately resolved, i.e., in roughly 
10(A/1.6/im)(D t ei/10m) _1 years, where A is the wavelength of the observations and D tc \ is the 
diameter of the telescope used, assuming the observations are diffraction limited (for a discussion 
of detecting light from a main sequence lens, see the Appendix). 

The results of the Bayesian analysis are shown in Figure [6l We find that if the lens is a main 
sequence star, its mass is Ml = 0.43^^ M and its distance is Dl = 7.15 ±0.75 kpc (median and 
68% confidence interval). Hence the planet mass is m p = 2.4^ 'g Mj up . In the close solution, the 
projected separation is sharply peaked at r± = sDlOe = 1.0±0.1AU. However, the wide solution, 
which is not strongly disfavored, gives an alternative r± = 3.4 ± 0.4 AU. If we assume a ~ r±, the 
planet would have a period of ~ 1.5 or ~ 8 years. 



8. Discussion 



8.1. Implications for Planet Formation Theory 



The lens in MOA-2011-BLG-293 consists of a super- Jupiter orbiting a probable M dwarf. The 
projected separation of the planet from the star is at most a few AU, making it difficult to form 
in situ if the host is indeed an M dwarf. Core accretion theory makes a general prediction tha t 



massive Jovian planets around M dwarfs should be rare (jLaughlin et al 



2004: 



Ida fc Linll2005l l. 



While gravitational instability can form large planets around M dwarfs (|Bossll2006l ). these typically 
form farther out, so if the planet formed by this mechanism, it would either be required to have 
migrated significantly or the projection effects must be severe. In the Appendix, we discuss how 
Adaptive Optics (AO) observations can confirm the microlensing measurement of the host mass or at 
least place upper limits on the host mass that would confine it to the M dwarf regime. Additionally, 
it should be noted that this planet joins a growing number of massive planets orbiti ng stars likely 



to be M dwarfs discovered b y microlensing (see Appendix) and radial velocity (see iBonfils et al 



20 111 for a summary and also lJohnson et al.ll2011 



8.2. The Ongoing Importance of Followup 

Even with high-cadence surveys, followup data remain important for the interpretation of 
individual events. In this case, the event was high magnification with a faint source, undetected 
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at baseline. This meant that the time period when the event was observable, to surveys and 
followup, was very brief. We have shown that followup data can vastly increase the signal and 
provide redundancy in the light curve coverage, which protects against weather. In Section 15. 2[ 
we demonstrated that with the survey data, even though the break in the light curve is missed 
because of bad weather, the same solution is recovered for MOA-2011-BLG-293, albeit at much 
lower significance than when all the data are included (A^ 2 ~ 500 compared to A% 2 ~ 5400). This 
is somewhat surprising since it is conceivable that because of the degeneracies possible for central 
caustic type events, the loss of this feature would leave the models relatively unconstrained or 
allow alternative solutions to the light curve. The followup data are beneficial in this case because 
in addition to adding to the signal-to-noise, they trace the sharp feature seen in the model light 
curves, increasing our confidence that the correct model has been found. 

It is likely that without the real-time discovery of this event in the MOA data and the subse- 
quent high magnification alert from //FUN, the planet would have remained undiscovered as of this 
writing. The planetary anomaly only becomes detectable when data from all three survey telescopes 
are combined, which requires systematically reducing, combining, and searching all of the survey 
data, preferably using a difference imaging reduction in order to detect events with faint sources 
like this one. Routine reduction of all survey data is planned for the current OGLE/MOA/Wise 
survey, but has not yet been fully implemented. 

The multi-band data taken by followup groups can also be important for interpreting mi- 
crolensing events. In order to determine the physical characteristics of the lens from the microlens 
variables, we used the (V — I) color of the source to estimate its angular size, thus providing a 
physical scale for the lens system. We also used the (V — I) color to inform our choice of limb- 
darkening coefficients for our model. Additionally, .ff-band data are important for comparison to 
AO observations that may be used to improve constraints on the lens mass as discussed for this 
event in the Appendix, //-band data are routinely taken as part of followup observations at CTIO 
but not part of the planned surveys. The OGLE survey regularly takes U-band data every few days 
when the weather is good, and the Wise survey is planning similar observations for future seasons. 
Provided the event timescale is long, this will result in several points taken when the source is 
substantially magnified. 

For the present case, the highly magnified part of the event is brief, and OGLE only obtained 
one V band point over peak, which was taken deliberately as part of followup observations. Because 
the U-band data were only taken as part of followup, we need to consider the effect of excluding 
these data in the context of a p ure survey detection . In principle, the (V — I) color could be 



estimated following the method in lGould et al.1 ( 2010bl ). which takes advantage of the difference in 



the OGLE and MOA bandpasses to estimate the (V — I) color from -Rmoa — ^ogle- For this event, 
the uncertainty of this measurement from the fits to survey-only data is (Tr-i = 0.01 leading to a 
(V — I) uncertainty of ay-i = 0.04. In this case, the precision is not much worse than what we 
found from the standard technique using CTIO V- and I-band photometry, so the lack of survey 
U-band data would not have a major effect. 
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Finally, although this event clearly shows that a planet is detected at A\ 2 ~ 500 without 
followup data, this is smaller than the A^ 2 of any published microlensing planet, underlining the 
fact that low-significance signals have not been systematically explored. Events like MOA-2011- 
BLG-293 that are robustly characterized with followup data but with weaker signals in the survey 
data can be used to probe lower A% 2 signals and inform our understanding of the limits of what 
is detectable. For example, if we analyze a large number of events with a range of signal strengths 
in the survey data, we could determine a A% 2 threshold for which a known signal, seen in the 
complete data set including followup, can no longer be distinguised from the noise. Additionally, 
the results for central caustic events will help us understand a A% 2 threshold below which the 
model degeneracies mean that the "correct" solution, as determined from the complete data set, 
can no longer be recovered. We might also require that these events can be well-characterized, i.e., 
that degenerate central caustic models can be sufficiently disentangled so that the mass ratio is 
well-constrained^l. Understanding these thresholds will be important for analyzing large samples 
of events to study the planets as a population rather than individual discoveries. By analyzing a 
large sample of events similar to MOA-2011-BLG-293, we can empirically determine appropriate 
Ax 2 thresholds or investigate other statistics for both detecting planets in all microlensing events 
and characterizing them in central caustic events. 
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A. Appendix: Possible Constraints from AO Observations 

Of the 13 previously published microlensing planets, two are very likely to be super-Jupiters 
orbiting M dwarfj^l. In both cases, high-resolution imaging from space or the ground was needed 
to complete these determinations. OGLE-2005-BLG-071 has a mass ratio, q = 7.4 ± 0.4 x 10 -3 



1 We will leave the exact definition of this phrase to future investigations but suggest that it might be along the 
lines of constraining the mass ratio to an order of magnitude at 2a. 

Two more could be marginally included in this category 



Bennett et aUbood i lDong et alj|2009bh . 
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(lUdalski et al.1 120051 : bong et all l2009al ). similar to that of MOA-2011-BLG-293 analyzed here. 
Dong et al.l (|2009al ) subsequently combined Hubble Space Telescope (HST) and light curve data 
to determine the host mass Mj, = 0.46 zfc 0.04 Mq, implying that this was an m p ~ 3.5A/j up 



planet orbiting an M dwarf. iBatista et al.l (|201ll ) found an even higher mass ratio for MOA-2009- 
BLG-387 of q = 0.0132 ± 0.003. Their marginal detection or upper limit on the lens flux from 
8m-class adaptive optics (AO) observations allowed them to place an upper limit M < 0.5 M & 
(90% confidence) on the host, with a median estimate of Ml = 0.19 Mq, and so m p = 2.6Mj up . 

Based on our Bayesian analysis (Section 17. 2p . it is likely that MOA-2011-BLG-293L is another 
case of a super-Jupiter orbiting an M dwarf. As we now show, 8m-class adaptive optics (AO) 
observations could clarify the nature of the system. The main uncertainty in an AO measurement 
of the lens flux is the uncertainty in the source flux, since the two objects are superposed unless 
many years have passed since the event. Because an alternate model for the baseline flux exists 
in which the flux is not constant (see Section [4 . 2 H . we cannot be certain it is possible to reliably 
measure the blended flux. Nevertheless, for the sake of argument, we are going to assume that 
the baseline flux is stable, so the uncertainty in fg is dominated by the statistical errors from the 
MCMC procedure. This assumption can be tested after OGLE-IV has collected several more years 
of baseline data on the event. If these data confirm that the baseline is stable, the calculations in 
this section are applicable. If not, the limits estimated here are overly optimistic. 

As mentioned in Section [U the statistical error in fs is 4.6%, which is due primarily to 
correlations of the source flux with other model parameters, rather than errors from fitting the 
light curve to an individual model. Thus, the /f-band flux (in units of the instrumental scale 
of CTIO ff-band images) is known to e ssentially the same precision. Using standard techniques 



(jjanczak et al.l 120101 ; IBatista et al.l l201ll ) , this flux scale can be aligned to the AO flux scale to a 
few percent. Hence, the AO source flux fn,s can be known to about 7%. This means that excess 
light due to the lens can be securely detected at the 3 a level, provided that AH = Hl — H$ < 1.7 
mag. This quantity can be related to the physical properties of the lens and source by 

AH = H L - H s = AM H + AA H + ADmod (Al) 



where AM H = M HjL - M HtS , AA H = A u ,l - A H ,s> and ADmod = b\og(D L /D s ). Now, in the 
regime we will be considering, it is very likely that AAu = 0, but in any case AAh < 0, since the 
lens is in front of the source. Hence, we can conservatively ignore this term. The last term is 



ADmod = —5 log! 1 + 



D s 6 



AU kM l J 



-0.3 



D s 0.45 M< 







8kpc M L 



(A2) 



where in the second step we have inserted Equation (|14p and kept only the first term of the Taylor 
expansion of the logarithm^. Hence, the lens will be detectable provided 



M h ,l < M H ,s + 2.0, 



(A3) 



This approximation assumes that (Ds/Dl - 1) < 1, where Dl is implicitly embedded in 0b- 
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where 2.0 = (AH — ADmod) assuming D$ = 8kpc and Ml = 0.45 Mq. Prom its color and 
magnitude, the source is a late G/early K Bulge dwarf, so Mu,s ~ 4.2. Hence, all dwarfs Mh,l < 6.2 
are detectable, which corresponds to > 0.43 M . 

Such a detection or upper limit would not be absolutely secure. Detected light could in principle 
be due to a companion to the lens or source, or an unrelated star in this crowded, low-latitude 
Bulge field. Additionally, an upper limit on the lens flux could be attributed to a remnant host 
rather than a late M dwarf. Nevertheless, the probabilities of these alternative interpretations can 
be quantified, and it is important to do so in order to estimate the frequency of very massive planets 
orbiting M dwarfs. 
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MOA-2011-BLG-293 



Peak 




5744 5745 5746 5747 5748 5749 5750 5747.3 5747.4 5747.5 5747.6 5747.7 

HJD' HJD' 



Fig. 1. — The light curve of MOA-2011-BLG-293. The left-hand panel shows a broad view of 
the light curve, while the right-hand panel highlights the peak of the event where the planetary 
perturbation occurs. Data from different observatories are represented by different colors, see 
legend. The black curve is the best-fit model with a close topology (s < 1). The times are given in 
HJD'=HJD-2450000. 



25 



22 



_l I I I I I L 



_] I I I I I I I I I I I L 



-1.0 -0.5 



0.0 0.5 

(V-I)ctio 



1.0 



1.5 



Fig. 2. — Color-Magnitude Diagram of the event in instrumental (uncalibrated) magnitudes. The 
source is shown as the solid black point; the errors in the source color and magnitude are smaller 
than the size of the point. The centroid of the Red Clump is the open square with an X through 
it. The small points show the stars in the field, restricted to stars within 60" of the source because 
there is strong differential reddening on larger scales. 
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Fig. 3. — Observed MOA (open squares) and OGLE (solid circles) fluxes at baseline. The fluxes 
have been scaled by the source flux and adjusted so that the baseline is approximately zero. The 
solid line shows the expected flux from the model. The data have been binned by 30 days (right 
panel) and semi-annually (left panel). Data taken when the event is significantly magnified (hashed 
region: 5710 < t(HJD') < 5790) have been excluded. Note that the MOA data show significant 
variation at a level comparable to the source flux. 
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Fig. 4. — Caustic structure and source trajectory of the best-fit model of MOA-2011-BLG-293 in 
the source plane. The circle shows the physical size of the source, and its position at the time of 
the caustic exit (HJD' ~ 5747.5). The x-axis is the star-planet axis, and the origin is at the center 
of magnification. The scale of the axes is in units of the Einstein radius. 
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Fig. 5. — Comparison of point-lens fits (left) and planet fits (right) for "survey-only" data, (top) 
and all data (bottom). In both cases, the planet fit is clearly better than the point lens fit, but the 
difference is more significant when followup data are included. Note that for "survey-only" data 
the OGLE data have been thinned out to reflect the typical survey cadence. 
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Fig. 6. — Probability densities for the lens (host) as function of its distance (top) and mass (bottom). 
The vertical scale is set so that the total area under each curve is equal to one. Masses Ml > 
1.2 Mq are excluded by the flux-alignment constraint on the lens brightness (bold vertical line). 
The 68% and 90% confidence intervals about the median are indicated by the shaded regions. 
The discontinuities in the slope of probability distribution for the lens distance arise from overlap 
between the Disk and Bulge stellar density distributions. From the Galactic model priors, there 
is a significant probability that the host is an M-dwarf. High-resolution imaging could confirm or 
contradict this by direct detection of the lens light (see Appendix). 
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Table 1. Data 







Error Renormalization Coefficients 




Observatory 


Filter 


k 


6min 


Ndata, 


OGLE 


I 


1.75 


0.01 


274 a 


MOA 


MOA-Red 


1.25 


0.0 


78 b 


CTIO 


I 


1.56 


0.0 


63 


Wise 


I 


1.57 


0.0 


49 


Weizmann 


I 


1.74 


0.0 


54 


CTIO c 


V 






9 



a A ? data after binning. 

b A r data after binning. Restricted to 5743.5 < t(HJD') < 5749.5. 

c These data were not used in light curve modeling. They were only used to determine the color 
of the source. 

Note. — The properties of each data set are given along with the error renormalization coefficients 
used to rescale the error bars (see Sec. 12. ip . 
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Table 2. Lightcurve "Invariants" 



Model 




/lim 




qtE 




(days) 




(days) 


(days) 


close 


0.0756(5) 


9.81(6) 


0.0355(3) 


0.115(2) 


close with 


0.0754(5) 


9.82(7) 


0.0355(3) 


0.115(2) 


flux constraints 










close with 


0.0748(7) 


9.95(11) 


0.0355(7) 


0.113(2) 


parallax 










wide 


0.0754(5) 


9.85(7) 


0.0355(3) 


0.116(2) 


wide with 


0.0754(5) 


9.85(7) 


0.0355(3) 


0.116(2) 


flux constraints 










survey only 


0.075(2) 


10.1(3) 


0.039(3) 


0.109(7) 


survey with 


0.076(2) 


9.9(3) 


0.040(3) 


0.110(8) 


flux constraints 











a /i im = /s,ogle/'"0) where /s,ogle = 1 corresponds to a magnitude I = 18, so /n m has units of 
flux in this system. 



Note. - Comparing the invariants of the lightcurve (Sec. I4.ip shows that they are robustly 
measured both in terms of their uncertainties and their variation among models. 
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Table 3. Model Parameters 



Model 


9 

X 


to — 5747. 
(HJD') 


Uq 


(days) 


P 


a 
(°) 


s 


9 




7TB, E 


fS,Wise 
fs.CTIO 


fs.Weizmann 
fs,CTIO 


close 


658.9377 


0.4935(7) 


0.0035(2) 


21.67(96) 


0.00164(7) 


221.3(5) 


0.548(6) 


0.0053(2) 


o.(.) 


o.(0 


0.979(9) 


1.09(2) 


close with 


662.0860 


0.4935(6) 


0.0035(2) 


21.75(95) 


0.00163(7) 


221.3(5) 


0.548(5) 


0.0053(2) 


o.(0 


o.(0 


0.990(4) 


1.08(1) 


flux constraints 


























close with 


655.5644 


0.4924(9) 


0.0035(2) 


21.24(95) 


0.00168(8) 


221.5(6) 


0.552(6) 


0.0054(2) 


1.7(1.1) 


-2.4(1.5) 


0.94(2) 


1.04(3) 


parallax 


























wide 


662.8497 


0.4931(7) 


0.0034(2) 


22.49(98) 


0.00158(7) 


221.1(5) 


1.83(2) 


0.0052(2) 


o.(.) 


o.(0 


0.98(1) 


1.08(2) 


wide with 


665.9169 


0.4931(6) 


0.0033(1) 


22.64(98) 


0.00157(7) 


221.1(5) 


1.83(2) 


0.0051(2) 


o.(0 


o.(0 


0.988(5) 


1.07(1) 


flux constraints 


























survey only 


497.3160 


0.492(1) 


0.0038(2) 


19.8(1.0) 


0.0020(2) 


218(1) 


0.55(1) 


0.0055(4) 


o.(0 


0.(.) 






survey only with 


498.8901 


0.493(1) 


0.0038(2) 


20.0(1.0) 


0.0020(2) 


218(1) 


0.55(2) 


0.0055(4) 


o.(0 


o.(0 






flux constraints 



























W Note. — The mean and root mean square errors for the parameters of each model are given along with the x 2 f° r that model. 

The fits with "survey only" use only a subset of data representative of what would have been obtained without additional 
followup. Note that the parameters of these fits are very similar to the parameters of the other fits, but with slight increases in 
their uncertainties. 



